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Abstract: We discuss the strong coupling of heavy mesons to a pion ^, in the 
heavy quark hmit. This quantity is quite remarkable since its values as estimated by 
different methods (various quark models and the QCD sum rules), are surprisingly 
different. The present quark models are mostly based on free spinors and their pre- 
dictions depend crucially on the choice of the light quark mass. We propose a quark 
model based on the Dirac equation in a central potential, which gives a more refined 
description of Dirac spinors. We show that within such a Dirac model, the value of 
'g is stable and large: ^ = 0.6(1), where we assume no quark current renormalization 
{{gA)q = !)• Such a large result is strongly constrained by requiring that the model 
parameters fit the spectrum; we show that this implies a large "effective" light mass. 
It is also supported phenomenologically by a similar situation with heavy baryons, 
as well as by experience with nucleon (if one invokes additivity). We also calculate 
the couplings to heavy meson excitations, and show that the Adler-Weisberger sum 
rule is well saturated by a few levels (in contrast to the case of small g). We dis- 
cuss uncertainties of our approach, and rise several questions which remain to be 
answered. The main mystery is the large, unusual discrepancy with QCD sum rules 
for ^, whereas a good agreement is found for orbital excitations. 
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1. Introduction 

A precise value of the coupling of pion to heavy-light mesons {gB*Bn, go'Dn, or 'g 
which is (in the heavy quark limit) related to Qb'Btt as, Qb'Btt — C^^mB/ fn) g), has 
been often needed in phenomenology during the last decade (see Ref. flj, for review). 
Qb'Btt is particularly important, since it can be helpful in the analysis of the semilep- 
tonic B ^ 7T decay. Namely, the exclusive modes, B 7r(p)£z/£, represent the most 
promising way to extract the value of the least known CKM matrix element, \Vub\- 
For that, the essential theoretical input are form factors, describing the hadronic 
matrix elements in the whole physical region. The calculations which incorporate 
QCD in their first principles (lattice QCD and QCD sum rules) cannot fully solve 
the problem, i.e. they do not cover the whole physical range of g^, in particular 
the very small recoil region, near q'^ax- simple complementary theoretical idea is 
the one of the nearest pole (a vector meson) dominance (VMD). According to the 
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extensive study in Ref. 0, in the heavy quark hmit, the form factor f^^^{q^) is 
actually dominated by B*, when close to ql^ax- Then, one has 



/ rriB* 1 q 



m 



B* 



where fs* gB*BTT^B* /"i, represents the residue of the form factor at that particular 
(nearest) pole^. It includes gB*Bn, whence a main interest in evaluating this quan- 
tity. Of course, the complete functional dependence for this form factor is different 
from VMD, i.e. for lower (larger recoils, (f), where the contributions of higher 



excitations become substantial, it becomes : 



^ , = i ( fB'3B^B.ms^ + V , (1.2) 

2 \ m^, — q^ ^ m^, — q^ J 



where /b* Qb^b-k ^ib^-, are the Svr couplings to the the higher excitations, -B*. For the 
complete treatment, one should also account for multiparticle states, which cannot 
be incorporated in our simple approach. 

There cannot be a direct experimental indication on Qb'Btt, because there is 
no phase space for the decay B* Btt. The available experimental results for 
D* Dn ^, can be related to gB*BTr, through the heavy quark symmetry {HQS). 
Still, the upper limit that can be obtained in this way is not very constraining (see 
Tab. [H). One must then turn to theory. 

A first lattice approach related to the above considerations would consist in 
comparing the lattice form factor data near zero recoil to VMD in the case of D — tt 
(which is directly accessible on the lattice), and estimate Qd'O-k and thus gB*Bn- In 
this way, using the APE data Q for fD{f=o)-^TT{pi t)(^^)' "^^ere piat is the smallest 
momentum that can be given to a particle on the lattice, one obtains ^ = 0.41(6). 
It is important to mention that the appropriate chiral extrapolation was performed. 
Results have also been obtained by working on B —y n lattice data, which are 
obtained by extrapolation in the heavy quark mass 0. 

Very recently, a direct preliminary determination of qb* b-k on the lattice has been 
attempted by UKQCD 0. The value they reported, g = 0.42 (4) (8), is compatible 
with the one extracted from the form factor close to zero recoil. 

We turn now to the other systematic theoretical methods for calculating the 
couplings: QCD sum rules (QCDSR) and quark models (QM). A new problem then 
arises: they lead to very different results for a seemingly simple quantity. One would 

2We take U = 131 MeV. 

3r(i:)*+)(e^p ) X B{D*+ r)^T^+){'^xp.) ^ qq^2) keV, amounts to the bound < 0.57. 
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like to understand why. The main resuh we present here, concerns the couphng ^. 
A short survey of the hterature indicates that the value of ^ is far from being estab- 
lished {see Tab. |I]). The overall conclusion drawn from different QCDSR calculations 
is that the coupling gB*Bn is small. In the heavy quark limit, their typical values are 
in the range ^ — 0.15 -h 0.35. On the other hand, the values proposed by different 
quark models are quite dispersed, ranging from ^ = 1 in nonrelativistic QM [p^] , 
down to ^= 1/3, in certain partially relativistic models (see Ref. [0, for instance). 
Such discrepancies are very unusual, especially for ground state properties. 

From the physical point of view, a question arises when one tries to reconcile 
the small values for ^with the Adler-Weisberger sum rule (AWSR). The experience 
with the coupling to vrA^, taught us that lowest level resonances give the essential 
contribution to the saturation of the sum rule. This certainly would not be the case 
for ttB scattering if the value of ^ were very small. Another general problem for a 
small ^ is the additivity property which is expected to relate the axial couplings of 
light quarks inside the various hadrons (see below for more discussion). Based on 
the study of the decay of charmed baryons [^], where the coupling is known exper- 
imentally, it can be inferred that, in general, there is no such a strong suppression 
(of?). 

We want to rediscuss the problem within the quark model approach, to see 
whether one can get more stable predictions, restore the saturation of AWSR as 
discussed above, and retrieve the additivity. This is achieved by working in the 
Dirac equation framework^. After presenting the basic definitions and motivations 
in Sec. H, we discuss the potential and sketch the calculation in Sec. |^. Discussion of 
the results are presented in Sec. ^. We end with a short conclusion in Sec. |[ 

2. Dirac equation, pion emission and Adler— Weisberger sum 
rule 

A nice feature of the model based on Dirac equation in a central potential is that it 
treats relativistically the quark internal velocities, with a certain amount of precision. 
In particular, it includes the Z graphs and appreciably modifies the estimate of the 
small components as compared to free spinors. The latter is crucial for the estimate of 
quantities such as axial couplings. It presents exact additivity for the current matrix 
elements. This model was initiated by Bogolioubov who used it to calculate qa 



- nucleon axial coupling (related to the vrA^ scattering). In the soft pion limit, one 
encounters the famous AWSR: 

9\- - K^+P) - a(7r-p)] = 1. (2.1) 



^While this paper was in writing, a similar strategy to calculate the pionic couplings was also 
considered in Ref. Roberts. However the subtleties of our models are quite different. 
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Lattice QCD 



Ref. 



0.42(4)(8) 



Direct method by studying 
the {B\A^{0)\B*). 



QCD Sum Rules 



Ref. 

Ref. i 
Ref. i 
Ref. |10| 



Ref. Ill 



0.39(16) 1 


S-R + soft pion limit; 


0.21(6)'^ J 


"■'> Included 0{as) corrections. 


0.15(3) 


Double moment SR. 


0.2 ^ 0.4 


SR in external axial Held. 


0.32(2) 


Light cone SR (calcul. of gB*Bn)- 


0.21(7) 


Light cone SR (calcul. ofgB'B'ir) 



Quark Models 



Ref. 112 1 
Ref. ig 
Ref. PI 



Ref. PI 
Ref. 1^1 



1.0 

0.60(5) 

0.57 
0.46(4) 

1/3 



Non-relativistic QM. 
Light front QM-I. 

Light front QM-U. 

Model which includes HQS 
and approximate xS. 

Partially relativistic model. 



This work 



0.6(1) 



With the preferred set of parame- 
ters: g = 0.61 (see the text). 



"Experiment" 



Ref. |T7| 



Ref. 118 



< 0.76 



0.27 tl 
0.76 tl tl 



Experimental upper 
bound (see footnote 1 ). 

From experimental data and 
accounting for the HQS and xS. 



Table 1: Results for g as obtained from the lattice QCD, various QCD sum rules and 
several quark models. In Ref. / [I^/ , the lagrangian approach (reviewed in j^) and the 
experimental results for 3(0*^^ D^j) and 3(0*^^ D^n'^) allowed the extraction 

of g, but two solutions were found. Along the same lines, but using the experimental 
B{D^ K^iv() and fos^ ^'^o distinct solutions for g were obtained also in Ref. / [I^ . 

Bogolioubov obtained the result: qa = 5/3 (1 — 26), where a positive 5 is a 
measure of the small components of Dirac spinors. With large quark velocities found 
in hadrons, the non-relativistic result = 5/3 is significantly reduced, and a good 
agreement with experiment (1.267(4) ||^) was reached. In the saturation procedure 
with the nucleon in the game, one should pay attention to quark's isospin and then 
sum over all quarks. In the case of meson with a very massive flavor, the situation 
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is even simpler since we deal with only one light quark (heavy quark decouple due 
to the heavy quark symmetry). Thus, the situation is very favorable for the use of 
the Dirac equation, which precisely corresponds to a particle in a static center of 
force. This is not only true for the description of the spectrum, but also for a quan- 
tum emission from the light quark (radiative transitions, pion emission/absorption)^. 

The exact validity of the AWSR in this (Dirac) formalism has already been 
demonstrated in Ref. We now recall some basic elements. Consider the axial 
current 

= ^(a^)7^75yg(x) , (2.2) 

and its corresponding charge 

Q\{n) = j d^x 6{n ■ x) n'^A;(x) , (2.3) 

where r* are ordinary Pauli matrices, whereas n = (n'',n) is a timelike unit vector. 
In the null-plane limit {\n\/n^ 1), the matrix element between one-quark states, 
X^^(n) = {a \Qa{^)\ P)y describes the pionic transition (m^r = 0): 

Kpin) = J d'x 5{n ■ x) ^o.{x){n ■ ^)^J-mp{x) , (2.4) 

where '^i{x) represent a complete set of solutions of the Dirac equation, which we 
take to be the eigenstates (\E'j(x) = \Efj(r)e~*'^'*), whereas n is chosen to be coUinear 
with the pion four-momentum 

K,{n) = J df ^r„(f) (^70 - 1^7) 75y *Mr) e^^^^"^")'^-''/'^" (2.5) 



The last step is to use [Q\{n) , Q^{n)~\ = 21^{n), demonstrated to be valid in the 
model in the null-plane limit (|n|/n° — > 1). Iz{n) is the corresponding vector 
charge. One finally arrives at 

Xq,5 X^^ = 5ap , (2.6) 



i. e. the Weinberg- like, matrix form of the AWSR . By specifying the null- vector 
n as = — >• 0, we have rewritten (Of) as: 



= j df^>l{r){^, - ^t)^>,{r) e-^^' , (2.7) 



^Thus for the processes when the recoil of the heavy quark can be neglected. 
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where S = 1 (8> <?, as usual, and with r matrices dropped out thanks to the additivity 
feature of the Dirac modeL As we aheady mentioned, a heavy-hght meson is (in the 
heavy quark hmit) the one hght quark system, and we can set \a) = \(3) = \B), and 
write 

J2\Xbs\' = 1. (2.8) 

6 

This is the form of the AWSR that we are interested in. X's will be directly related 
to the pionic decay widths through the Weinberg's formula (m^r = 0): 

r(J^^7r) = ^— VlXj^^p, (2.9) 

with sums over initial and final polarisations and over final charge states, and an 
additional factor 1/2 for neutral pion (/^r = 131 MeV). This will be the main 
formula for our investigation. 

We must emphasize at this stage, that a basic assumption we make is that the 
QM current is identical to the bare one (in terms of quark operators). Namely, 
the QCD radiative corrections in principle introduce a {gA)q factor for quarks P^], 
representing the phenomenological renormalization for constituent quarks. We do 
not introduce such a factor for two reasons. First, we want to remain within the 
point of view of saturating the AWSR by states composed of valence quarks only. It 



is rather obvious that the sum (|2.8| ) will no more equate 1, if the axial quark current 
is multiplied by {gA)q 7^ 1- On the other hand, if one includes radiative corrections, 
then one should also introduce the states with gluons, to keep the consistency with 
AWSR. Second, such a factor {gAjq notably different from 1, does not seem to be 
required to explain the (similar) charmed baryon decays, or the nucleon {gA)N (see 
Sec. [4.31 ). At least, it could not be very small. Nevertheless, a reduction factor, not 
too different from 1, is by no means excluded. 

In phenomenological Lagrangians (see [|l], ^ and references therein), the pionic 
transitions are basically described in terms of invariant coupling constants. This has 
several advantages: the phenomenological couplings incorporate explicitly Lorentz 
invariance and the consequences of the properties for the decay; they may also 
include the consequences of the HQS and chiral symmetry {xS). From the HQS, we 
know that the wave function of a heavy meson is independent of the heavy flavor and 
its spin. Consequently, the total angular momentum of the light degrees of freedom 
( i.e. of a light quark) ji, is a good quantum number and heavy mesons come in 
degenerate doublets with spin J = ± |. The doublets of interest here, are: 

5(iO-) , 5*(a-) ] ; [ i?*(iO+) 5;(il+) ] ; [ B^i|l+) , i?;(|2+)] (2.10) 

with the notation, j^J^. From the parity conservation, one can easily see that: 
i) the transition from one meson to the other within a doublet, occurs with the 
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emission (absorption) of the P-wave pion; ii) between B and Bq, as well as between 
B* and Bi{B[), the transition occurs with the S-wave pion; iii) between two mesons 
of next to neighbouring doublets, with the D-wave pion (when allowed), and so 
on. The covariance properties of these transitions are parametrized by the invariant 
couplings, QBiB-K) with B denoting here generically heavy-light meson, and Bi stands 
for B* , Bq, i?2 and (in general) for their radial excitations too. We will henceforth 
leave aside the l"*" states to avoid the complications of mixing, which anyway cannot 
be described by the Dirac equation. 

{Btt I B*) = {q^e^')gB*B. 
{B n \ Bq) = gB*BiT 




{B n\B;) = ^ - {e,,q^qn Qb^b. (2.11) 

is the pion four momentum, whereas e^{e^j_y) is the polarization vector (tensor) of 
5*(i?2)- In this definition gB*Bn = gs^oB-w+y is related to the other charge combina- 
tions by isospin: 

gB*Bn = gB*OB-n+ = -fi's—fiO^- = V^gf^.O^O^O = -\/2gB*+B+nO ■ (2.12) 

gB^B-K and gB^Bn have the mass dimension [m] and [m~^], respectively. Obviously, 
the above definitions apply in the same way to D and B mesons. The corresponding 
decay widths are: 



|3 



o A ^1 ID Bit 2 ' 

247r m^. 



^0 



V{Bf B^n-) = gl.^^ ^ , (2.13) 

407r 2 jYi^ 



By writing this generically, we have 
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\2l+l 



V{Bf fi+vr-) = gl,^^ (2.14) 

STT [ijB* + 1) ^B* 

* i 

with I = Jb* being the spin of the decaying meson. If one includes the constraints 
of HQS and xS in the phenomenological lagrangians the pionic transitions 

are described (at lowest order) by the new, dimensionless couplings: {g = g), h, h', 
which are the crucial parameters of the approach and are independent of the heavy 



^Remark again that unlike for r{D* — > Dtt), the case, B* Bn, is kinematically forbidden. 
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mass. In fact, the X's which we defined above (|2.7|) , are essentially these couplings^. 
To obtain the exact relations with our X's, we compare the corresponding widths in 
the heavy quark limit. In terms of the new couplings, the expressions ( p.l3|) read ^ 



^2 



Gicf^ TUB' 

r(fi*° B^n-) = J^'l ^ I g 1^ . (2.15) 

On the other hand, we know that the soft pion emission expressed in terms of 
X's is (PI) 

rW°-a^O = ^^J^|x.^ (2.16) 

with Xj calculated for the u quark and zero polarisations, so that we can identify 
(up to phases): 



2 a/ vfiB niB* ^ ^ 

9b*b-k = 7 9 ; 9 = 9 = Xb* , 

I mB ,2 2 N ^ > V 

9b*,b. = ^ ^ (m^. -^b) J^-, h = Xb* , 

[Y ^ niB* rriB h' h' _ [Y X^j 

9BiB.-^Y U A^' T^-\IYAE- ^^-^^^ 

The kinematical factors in the first column are kept to their physical values, to 
(somewhat) account for finite mass effects, although the higher (in l/mq) couplings 
are being neglected. On the other hand, the second column is obtained consistently 
in the heavy quark limit, since the couplings g, h, h' are the only dominant in this 
limit. X's are to be obtained from the Dirac equation (in the itlq oo limit), where 
we will also assume m^r = 0, i.e. uj^^ = \ q \ = AE [AE is the difference between the 
final and initial Dirac energy level). ~ 1 GeV, is the chiral symmetry breaking 
scale. Admittedly, the choice we make is somewhat arbitrary, although conformal 
to our experience with quark models. Another reasonable procedure which can be 
used, consists in fixing the matrix element cu^ X/Z^r, to its limiting value, AE X//^, 

^Indeed, the above X's are defined for the Dirac equation, and will be valid in the heavy mass 
limit (mg — > oo) of the theory. 

*For the B* — > B transition, our formula derived exactly from the Lagrangian differs by the 
factor mB/mB* from the usual expression. 
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and to take only the phase space factor at its physical value. The resulting widths 
are somewhat different which indicates the uncertainty in the results. 

From the above relations, we see that (in the rriQ oo limit) the AWSR (|2.^ ) 
can be rewritten in the following form: 

ig = gr+h' + -im'^ + -- - + '^'^'^^ + ^"^^^■^'^^^^■^^^ =1.(2.18) 
3 excitations states 

" . ' 

orbital excitations 



3. Calculations with a specific class of Dirac equation 

3.1 Potential and solution of Dirac equation 

To make use of (|2.5|) , we have to find eigenstates of the Dirac equation: 

(^a-p + m(3 + V{f) - E^^{f) = (3.1) 

and this is the part where we have to fix the model. Models based on Dirac equation 
differ in the form of the potential V{f). In particular, the differences come from the 
Lorentz structure of the confining part of the potential. Here, we will take the linear 
confining potential Lorentz scalar. 

V(r) = -— + (3 (ar + c) . (3.2) 
r 

The Coulomb constant k, the string tension a, and the constant c (which is included 
in the scalar-confining part ^), are the parameters of the potential. Note that only 
the combination (m + c) enters the equation ( p.l|) . Therefore, there will be only 
three parameters in the model: a, k and (m + c). It must be emphasized that 
neither m, nor (m + c) corresponds to the constituent quark mass, so that neither 
should be constrained to remain around 300 MeV; (m + c) may even be negative (see 
helow). In Ref. this model was discussed with various forms of the potential. On 
theoretical grounds (see P^, PS[), the Lorentz structure of the confining potential is 



a delicate issue. There is no reason to take it as a vector, since it cannot correspond 
to the exchange of one single gluon, even if dressed by very large self-energy effects. 
Indeed, general considerations for the heavy quark-antiquark potential lead to a 



spin-orbit force with a sign opposite to the one predicted by a vector exchange [2£ 



Yet, this does not allow to say that the confining potential is the Lorentz scalar. 



^If a constant c were included in the vector part, the resulting effect would be an overall shift of 
the spectrum. 
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The scalar is only the simplest choice which is in agreement with the quarkonium 
phenomenology, and in particular with the spin structure (which is very sensitive 
to the Lorentz assignment of the potential). Therefore, it is reasonable to calculate 
our X's adopting such a potential (|3.2| ). If the constant c is included in the scalar 
confining part (and there is a priori no reason to discard it), only the combination 
(m + c) (and not only m) enters the game. In other words, we cannot distinguish 
between the 'quark mass' and the constant in the potential. This is not worrisome, 
because m does not play the role of the constituent mass, in the framework of the 
Dirac model (see below). 

For the solution of the Dirac equation, we will heavily rely on the work of the 
Wisconsin group on the D, Dg, B, Bg (spin averaged) states where they also 
discussed the way how to numerically treat the Dirac equation with the potential. 
We differ however on the questions such as the string tension, the quark mass and 
the constant in the potential (which we allow to be non-zero). 

The general solution of the spherically symmetric problem at hand, can be writ- 
ten as: 



if) 



3m\ 





II 




^3^7™ (^) ) 







(3.3) 



yimi^) denote usual spherical harmonics, where k stands for the Dirac quantum 
number {k = I for I = j + and k = —{I + 1) for I = j — |). After inserting this 
into the Dirac equation, one arrives at the system of differential equations. 



fm + ar + c-^-E 
k 



\ 



+ 



A 

dr 



k 

r 



Hr^ r 



\ 



— (m + ar + c) 



K 

r 



E 



9jir) 



0. (3.4) 



N 



As it was proposed in , we look for the solutions of the radial part 

N 
i=l 

in the pseudocoulombic basis of functions: 



- cj^^ 'Pirn 

i=l 



(3.5) 



8(i- 1)! 
(i + 2/ + l)! 



/5i 



(2/3r)' e"^'- Lf+^{2(5r) 



(3.6) 



For a certain A^, one diagonalizes the hamiltonian and gets A^ positive and A^ negative 
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energy solutions. The eigenvalues Ei, should not depend on the variation parameter 
(3, which is actually the case if is large enough (in our case it turns out to be 
for N > 15). In order to fix the parameters of the potential, one has to compare 
the predicted spectrum with the experiment, which is not a trivial task. The Dirac 
equation has the properties of heavy quark mass limit, i.e. the mass of a heavy-light 
meson is mn = uiq + Ei + 0{l/mQ). If we neglect the terms (9(l/mQ), rriQ can be 
adjusted trivially for each quark species to fit the ground state, or more precisely, 
the spin-averaged (0^, 1^) ground state. However, this approximation appears to be 
too crude, as we show below. 

Let us reconsider the problem of comparing the experimental spectrum with the 
predictions of the Dirac equation. Note that in the latter only the spin-orbit force 
(relative to the light quark) survives, whereas the spin-spin force is not present^°. 
What is really predictable and comparable with experiment is: AE^, the mass 
difference between the spin average of the {L = l,j£ = 3/2) states and the spin- 
averaged (0~, 1") ground state (extrapolated to the heavy quark limit), and AEg o., 
the spin-orbit splitting between (L = = 3/2) and (L = = 1/2) states (in 
principle also to be extrapolated). 

(A£'^)^''^P) is known, modulo spin averages and an extrapolation in l/mq. By 
using D and B states {i.e. AEl and AE\) and extrapolating as 

AE^ = AE^ + (3.7) 

one finds that the coefficient A is sizable A ~ 0.16 GeV^. That means that for D 
mesons, the corrections are as large as ~ 25%. 

As for the L = = 1/2 states, they are yet to be detecteded, although an 
interesting constraint on the spin-orbit splitting can be deduced from existing data. 
Namely, the observed (broad) B** bump cannot be a priori ascribed to the (narrow) 
je = 3/2 states only. That means that the L = l,je = 1/2 states are also present 
in B**, and their widths should be around 100 MeV (which may be the explanation 
why they were not isolated). This fact gives an important, additional constraint on 
the spectrum: the splitting between ji = 3/2 and je = 1/2 states must be small, i.e. 
such that only one bump stays 'visible'. 

Thus, the two firm experimental facts which we can compare with the predictions 
of the Dirac equation (and thus fix the parameters of our model) are: the value of 
AE^, and the (small) spin-orbit splitting, AEg o.- 



The set of parameters which reproduced the spectrum in was: 



m 



300 MeV , a = 0.308 GeV^ , k = 0.579 , (3.^ 



-'^'^This implies that e.g., the L = 1 states with same are degenerate. 

"OPAL Coll. H, and DELPHI Coll. |§ reported T{B**) ^ 116(24) MeV and r{B**) 
145(28) MeV, respectively. 
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with constant c set to nought. The manifest problem of such a model is the large 
value of the string tension, which is normally expected to be, 0.15 ^ a < 0.20 GeV^. 
After inspection, we realized that the problem appeared because they fit the spectrum 
of the D states {^El = 0.47 MeV), while the Dirac equation is valid in the infinite 
mass limit. From the measured B** (L=l) states, we can infer AE^ ~ 0.38 MeV, 
which means that it should be even smaller when mg —>■ oo. Linear extrapolation 
in Eq. ( pI7| ) amounts to AE^ f« 0.35 MeV, close to the existing estimate In 
order to fit this {smaller) value of AE^, a much lower value of the string tension 
a is needed than that found by authors of Ref. ||2^. In other words, the string 
tension is very sensitive to the value of AE^. Since AE'^ < AE^, we can fit the 
spectrum with standard values for a. Reasonable values of the other parameters are 
—0.2 ^ (m- + c) < 0.3 and k > 0.6, the simultaneous choice of the three parameters 
must be correlated. Our preferred choice, which was used to get the results listed in 
Tabs. is: 



-{m + c) = 200 MeV , a = 0.20 GeV^ , k = 0.65 . 



(3.9) 



Parameters of the potential (Eq. |3^ ) 
a = 0.20 GeV^ 
K = 0.65 
-(m + c) = 200 MeV 



Dirac spectrum 

AE^ = 341 MeV 
AEs.o. = 1.3 MeV 
AE' = 494 MeV 



Table 2: Our preferred set of parameters used for the spectrum for which we tabulate: 
AE^ = E£)irac{L = l,k = —2) — E£)irac{L = 0) ; the small energy splitting within 
L = 1 states due to the spin-orbit forces in the mg = oo limit, is AE'^.o. = Eoiracik = 
—2) — E£)irac{k = +1), and AE' is the distance between the ground state and the first 
radial excitation in the same limit, also given by the Dirac equation. 



3.2 Expressions for pionic transitions 

Now we go back to our problem and calculate the pionic transitions from the obtained 
wave functions, solution of the above Dirac equation ( p.l| ). To end that, we take the 
pion emitted along the 2;-axis, and use the plane wave expansion in (PT7|), to end up 
with 



1=0 



21 + 1 



rff v^yf)(7^ - E,)^BAf^ji{qr)Yi'{r) (3.10) 
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In terms of the eigenf unctions of the Dirac equation, the couplings (for transitions 
from excited to the ground state) read as follows: 



X 



BB* 



|f{-l)|2_ i|o(-l)|2 

\Ji/2 I r^\yi/2 I 



X 



BBr, 



(1) _ 

Ji/2 yi/2 yi/2 Ji/2 



3o{qr) + fi/2^*fyl + 91/2 9[/i^ jiiqr) | r^dr 



V2 



.{-!)* .(-2) _ i^(-l)* (-2) 
Jl/2 J3/2 e,yi/2 93/2 



Jiiqr)- f[/2*9i/2 - 9[/2*fi/2^ hiqr) 



7 9[/2*9i/2h{qr) W'^dr. 



(3.11) 



For eigenvectors, the dependence on r is imphcit (/•^'' = fj^'{r) and g)''' = gf'{r)) 



' 3l 



in the above formulae. In a given approximation, the same expressions apply for the 
radial excitations too, where the corresponding radial functions (/ and g) are eigen- 
vectors of the corresponding eigenvalue, E'-^i2 = 795 MeV, of the hamiltonian (|3.4| ), 
in the basis (|3.5|) with = 30. Note again that we take q = \ q \ = AE'; there- 
fore, I (f I = 0, for D* Dtt transition. The resulting values for the X's and the 
couplings g, h, h', are listed in Tab. ^ whereas the corresponding decay widths are 
tabulated in Tab. El 



4. Discussion of the results 

4.1 Discussion of the numerical estimate of ^, comparison with QCD sum 
rules and experiment 

For I g* I = 0, the expression for X= ^ is particularly simple: 



X 



BB* 



in 



(-1)12 

/2 I 



with the normalisation condition. 



(4.1) 



I . (-1)12 
1/1/2 I 



+ \9[-;2^\' \ r'dr. 



(4.2) 



In the non-relativistic limit, the small component, g^^^ is zero, and therefore: ^ = 1. 
This result would correspond to the famous 5/3 for the nucleon. As we already 



13 



Ground state (j^L =1/2 0) elastic transition 



Preferred parameters 
a = 0.20 GeV^ k = 0.65, (m + c) = -0.2 GeV 



X 



1/20^1/20 



g = 0.607 



Another set of parameters with satisfactory spectrum 
a = 0.15 GeV^ k = 0.62, (m + c) = 0.3 GeV 



X 



1/20^1/20 



g = 0.728 



Transitions to (L = 1) excited states 

Preferred set of parameters 

l-^i/20-i/2il = 1^1 = 0.536 
^./20-3/2il = 0.234 , \h'/A^\ = 0.841 



Table 3: Transition matrix elements X's, and couplings of heavy mesons with a single 
pion. The basis of functions {31,3^) is used with N=30, and {gA)q = 1, as discussed in 
the text. 



explained, with the small components different from zero, the value for ^ is consider- 
ably reduced. In abstracto, it could be very small, but we obtain (with our preferred 
parameters) , 



^ ~ 0.61 , 



which is by far larger than the values advocated by QCD sum rules ^ ~ 0.2 0.3. 
We stress that this moderate suppression in our model remains such, even when 
we vary the parameters in their reasonable ranges (see Sec. We obtain, ^ ~ 
0.6 ± 0.1. This value corresponds to what is needed to get gA/dv < 5/3, i.e. 
closer to experiment |2l|. To better monitor the consistency, we also observe that 
if we adopt the free Dirac spinors and m = 0, we get a result of g = 1/3, 
but the success of Bogolioubov's result for gA/gv would be totally lost. In fact, 
in this argument, we invoke the general principle of additivity for the coupling of 
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gB*B7r = 49.1 


go'D-K — 18 


T{B* Btt) = NO 


V{D* Dtx) = 93.7 keV 


gs^^B. ^ 19.5 GeV 


gD*Dn ^ 8.2 GeV 


T{B* Bn) ~ 290 MeV 


T{D* Dn) ~ 320 MeV 


gs^Bn = 58 GeV-^ 


go^^Dn = 22.5 GeV-i 


T{B* S(*)7r) = 27 MeV 


T{D* D^*h) = 47 MeV 



Table 4: Decay widths, calculated in the way described in Sec. |^ ; we use the masses 
niD* = 2.4 GeV, mD* = 2.46 GeV, tub** = 5.732 GeV. Note that the widths are sensitive 
to the assumptions concerning excitation masses. 

light quarks within different hadrons, which is exact in the framework of the Dirac 
equation : the light quarks are considered to be independent and bound to a common 
center of force. In the nucleon, it is possible that the equal light quark masssituation 
implies smaller internal velocities than in the D mesons, which could invalidate the 
use of our Dirac model and the additivity. However, the additivity argument works 
well for light quarks in charmed baryons which are more analogous to D mesons, 
and one finds indeed that the suppression with respect to the naive model seems 
experimentally moderate [pO|]^^. 

As for the sensitivity of our (relatively large) result on the choice of parameters, 
we emphasize that we have quite clear conclusion that ^ > 0.5, regardless of the 
choice of parameters provided they are compatible with the spectrum, 'g decreases 
with (m + c), but the latter should be kept > —0.2 GeV to ensure that the spectrum 
is well reproduced. We repeat, in the framework of the quark model with Dirac 
equation, low values of ^, like 0.2 or 0.3, are completely excluded. 

On the other side, to keep ^ compatible with the experimental upper bound 
(see Tab. [ip, we have to ensure that (m + c) < 300 MeV^^. With (m + c) small 
and negative, the spectrum can be very well described, and we are safely below the 
experimental bound. This determined our preferred set (|3.9| ). We must note that 
0{l/mc) enhancement with respect to the rriQ oo prediction is expected. In [0], 

^^We stress again that gA/gv and g, should be also suppressed by the QCD radiative correc- 
tions [ p4| , and therefore the small component correction must not be too large. 

^^With (m + c) ~ 300 MeV, we have g ~ 0.73, and our numbers are compatible with the bound, 
but quite close to it. 
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it is shown that 




where Coo is the rriQ — > oo hmit of the l.h.s. Since fs'^/mq, deviates from its hmit 
by a negative Oil/mo) term, one sees that gB*B-n/^Q deviates from its hmit 'g/ f-^ by 
a positive 0{l/mQ) term. This is not in contradiction with our way of extrapolat- 
ing gB'Bn/TTT'Q (sQ. |2.17| ), siucc y/mBmB'/mq may be expected to be larger than one. 



4.2 Interpretation of the results of the Dirac equation model in compari- 
son with previous quark model calculations 

An intriguing question for quark modelists is why our findings are so different from 



the ones of Ref. |jT6| {i.e. 'g = 1/3). We already showed that this number can be 
obtained with a free Dirac spinor structure and m = 0, independently of the spatial 
wave function. This would correspond to 

/oo 
I I' ^"d^ ■ (4-4) 

However, the result is quite sensitive to the quark mass, m. On the contrary, in our 
case, the magnitude of the small components in the Dirac scalar confining potential 
is determined by an average effective mass, meff = {m + ar + c) (for simplicity, 
we disregard now the Coulomb potential) In other words, in our approach, the 
average {m + ar + c), and not just m or (m + c), corresponds to the idea of the 
constituent mass. This average is large for any reasonable choice of parameters fitting 
the spectrum ; e.g. for the lowest admissible (m + c) = —0.2 GeV and a ^ 0.2 GeV^, 
K e (0.6^0.7): 

a (r) ~ 0.2 GeV^ x 2.45 GeV-\ (4.5) 

which gives (m + ar + c) ~ 0.3 GeV, therefore, already a rather large number^^. 
Therefore, the integral over g^^^ (< 0.3), remains always much smaller than over 

fy2\ fo^^ the values of (m + c) compatible with the spectrum. In other words 
^ remains large. We reiterate that a smaller ^ is allowed only when (m + c) gets 
much below zero, which turns out to fail in the spectrum description. Therefore, the 
small values for ^, like 0.5 or smaller, are prohibited. 



^^We owe this remark to J.-C. Raynal. 

"'^^Larger values are obtained at larger (m + c), although (r) is smaller: e.g. for (m + c) = 0, we 
have TOe// = 0.42 GeV, and for (m + c) = 0.2 GeV, me// = 0.55 GeV. 
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4.3 Discussion of other picnic transitions, additivity and saturation of the 
Adler— Weisberger sum rule 

- Orbital excitations: First, let us notice that our infinite mass limit coupling h, is 
in a quite remarkable agreement with the findings of the QCDSR 0. On the 
other hand the value of h'/A^ is compatible with sum rules but somewhat larger. 
As for comparison with experiment : 1) As inferred from Tab. ^, the D** and B** 
{ji = 3/2) 2"*" widths, are reasonably well described by our model. Though, our 
is somewhat too large, which may be explained by large l/rriQ corrections in the 
D sector (see in this respect |^). 2) The predicted large width of order 300 MeV 
for Dq may explain why it is not seen. The experimental total width of the B** 
bump (~ 140 MeV) is smaller than our prediction, ~ 250 MeV, for the Bq ; this 
remains to be understood. One interesting feature of this (emission) model is the 
relative narrowness of the S wave (0^) decays, which were predicted to be much 
larger (around 1 GeV) by previous quark pair creation models |3^. This is actually 



due to the fact that our model contains an additional factor of uj^^ in the amplitude. 

- Recently, DELPHI claimed to have detected the very narrow, would be (2637), 
radial excitation |3^. Subsequent search by CLEO p9[^^ did not find the trace of 



such a 'peak'. In Ref. [^], the radial excitation interpretation was excluded. Our 



model is not well suited to discuss radial excitation decays because of the sensitivity 
to the nodes of the wave function EM. With this remark in mind, we give the result 



|X^./| = 0.237, ^ T{D*' ^ Dn) ~ 70 MeV . (4.6) 
In Ref. PDl, it was concluded that the experimentally measured E* — Ac + 



TT width is in agreement with a naive non-relativistic QM estimate, with a small 
reduction factor. The situation is similar to the one of the D meson, and one expects 
the validity of the additivity principle: the coupling of light quarks to pion inside 
is the same as in the D. This is certainly true for the model with Dirac equation. 
This observation makes us confident that the value of ^ is indeed large. 

- For the end, we verify the saturation of the AWSR. We rewrite ( |2.18| ) without 
radial excitations and multiparticle states: 

2 h'^ 

f + + -{AEf— + ... = 0.368 + 0.287 + 0.055 + ... , (4.7) 
3 

which means that the first pole and orbital excitations saturate more than 70% the 
AWSR for Btt scattering, with the B* contributing the most. The saturation seems 
to work with a small number of resonance levels (like in the case of nucleon), in 
contrast with what would happen if we took the result of QCD sum rules. 



^^N.B.: In this reference, presented are also the first prehminary results for a broad = 1 
charmed state: mo, = 2.461^;^^ ± .010 ± .032 GeV, T = 290t\°g^ ± 26 ± 36 MeV. 
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5. Conclusion 



Using a Dirac model with a scalar confining potential (which per se includes 
the property of additivity of pion couplings), we have obtained a prediction for 
^ = 0.6(1), which is well defined and stable under the variation of the model 
parameters (chosen in a way compatible with the observed hadron spectrum). 
This is somewhat different from the previous quark model calculations, which 
were very sensitive to a (not well constrained) light quark mass parameter. Our 
model also provides a good description of the decays of orbital excitations, and 
(presumably) of the charmed baryon decay, S* — > Ac + vr; it gives the right 



9a/ gv, for nucleon It saturates well the Adler-Weisberger sum rule with 
a few low-lying resonances, as observed for the nucleon. The result for 'g and 
consequently for gD*Dn and gB*Bn is definitely large: 'g > 0.5. 

Of course, the main concern is the complete discrepancy with QCDSR pre- 
dictions for ^. It is the first time that one observes such a large discrepancy 
between the two methods in a situation where both are expected by their pro- 
ponents, to work reasonably well. 

Our assumption of a scalar potential (which lacks a QCD basis), the use of an 
elementary quantum emission model for strong decay to a hadron, and finally 
the use of the bare quark current for an emission from a light quark (in spite 
of possibly large renormalization effects), are possible defaults of our approach. 
Hence, we cannot claim to have compelling and accurate predictions in a quark 
model. For instance, it is quite possible that one should have to introduce a 
reduction factor for the quark current {gA)q- Still, with all these elements in 
mind, it would be very surprising if our predictions turn to be totally away 
from the correct number (as the QCDSR results suggest). 

Since one can always doubt results of a particular model, we insist on our 
strongest argument: it would not be easy to escape our (large 'g) conclusion 
within any quark model, and to explain the magnitude of the known pionic 
couplings of ground state baryons, at the same time. Indeed, if we lessen the 
strength of ^, then a similar suppression in the pionic couplings in other ground 
state hadrons (due to approximate additivity for light quarks pion couplings, 
expected at least to work well for charmed hadrons) would be difficult to escape. 
Now, such a suppression is not actually observed. It is neither possible to 
introduce a small {gA)q in order to reconcile the quark model with the small 
result of sum rules, because it will be still more unavoidable to end with such 
a suppression in other pionic couplings, which is not observed. 

On the side of QCDSR, certainly the convergence of the various calculations 
towards a small value of 'g is very impressive. However, one should: 



18 



a) have a conceptual understanding of why the result is so small, as compared 
to the naive non-relativistic QM (~ 1/5 suppression factor!), while the 
latter works reasonably for the analogous Qp-ktt, and gx^Kn couplings {e.g. 
recent result of Ref. p]); 

b) explain why the heavy baryon widths apparently do not exhibit any par- 
ticularly strong suppression with respect to a naive non-relativistic pre- 
diction; 

c) understand how the AWSR is saturated, with the ground state contribu- 
tion being so small. 
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